\(\int \sin ^5(c+d x) (a+b \tan (c+d x)) \, dx\) [11]

   Optimal result
   Rubi [A] (verified)
   Mathematica [A] (verified)
   Maple [A] (verified)
   Fricas [A] (verification not implemented)
   Sympy [F]
   Maxima [A] (verification not implemented)
   Giac [B] (verification not implemented)
   Mupad [B] (verification not implemented)

Optimal result

Integrand size = 19, antiderivative size = 101 \[ \int \sin ^5(c+d x) (a+b \tan (c+d x)) \, dx=\frac {b \text {arctanh}(\sin (c+d x))}{d}-\frac {a \cos (c+d x)}{d}+\frac {2 a \cos ^3(c+d x)}{3 d}-\frac {a \cos ^5(c+d x)}{5 d}-\frac {b \sin (c+d x)}{d}-\frac {b \sin ^3(c+d x)}{3 d}-\frac {b \sin ^5(c+d x)}{5 d} \]

[Out]

b*arctanh(sin(d*x+c))/d-a*cos(d*x+c)/d+2/3*a*cos(d*x+c)^3/d-1/5*a*cos(d*x+c)^5/d-b*sin(d*x+c)/d-1/3*b*sin(d*x+
c)^3/d-1/5*b*sin(d*x+c)^5/d

Rubi [A] (verified)

Time = 0.10 (sec) , antiderivative size = 101, normalized size of antiderivative = 1.00, number of steps used = 8, number of rules used = 5, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 0.263, Rules used = {3598, 2713, 2672, 308, 212} \[ \int \sin ^5(c+d x) (a+b \tan (c+d x)) \, dx=-\frac {a \cos ^5(c+d x)}{5 d}+\frac {2 a \cos ^3(c+d x)}{3 d}-\frac {a \cos (c+d x)}{d}+\frac {b \text {arctanh}(\sin (c+d x))}{d}-\frac {b \sin ^5(c+d x)}{5 d}-\frac {b \sin ^3(c+d x)}{3 d}-\frac {b \sin (c+d x)}{d} \]

[In]

Int[Sin[c + d*x]^5*(a + b*Tan[c + d*x]),x]

[Out]

(b*ArcTanh[Sin[c + d*x]])/d - (a*Cos[c + d*x])/d + (2*a*Cos[c + d*x]^3)/(3*d) - (a*Cos[c + d*x]^5)/(5*d) - (b*
Sin[c + d*x])/d - (b*Sin[c + d*x]^3)/(3*d) - (b*Sin[c + d*x]^5)/(5*d)

Rule 212

Int[((a_) + (b_.)*(x_)^2)^(-1), x_Symbol] :> Simp[(1/(Rt[a, 2]*Rt[-b, 2]))*ArcTanh[Rt[-b, 2]*(x/Rt[a, 2])], x]
 /; FreeQ[{a, b}, x] && NegQ[a/b] && (GtQ[a, 0] || LtQ[b, 0])

Rule 308

Int[(x_)^(m_)/((a_) + (b_.)*(x_)^(n_)), x_Symbol] :> Int[PolynomialDivide[x^m, a + b*x^n, x], x] /; FreeQ[{a,
b}, x] && IGtQ[m, 0] && IGtQ[n, 0] && GtQ[m, 2*n - 1]

Rule 2672

Int[((a_.)*sin[(e_.) + (f_.)*(x_)])^(m_.)*tan[(e_.) + (f_.)*(x_)]^(n_.), x_Symbol] :> With[{ff = FreeFactors[S
in[e + f*x], x]}, Dist[ff/f, Subst[Int[(ff*x)^(m + n)/(a^2 - ff^2*x^2)^((n + 1)/2), x], x, a*(Sin[e + f*x]/ff)
], x]] /; FreeQ[{a, e, f, m}, x] && IntegerQ[(n + 1)/2]

Rule 2713

Int[sin[(c_.) + (d_.)*(x_)]^(n_), x_Symbol] :> Dist[-d^(-1), Subst[Int[Expand[(1 - x^2)^((n - 1)/2), x], x], x
, Cos[c + d*x]], x] /; FreeQ[{c, d}, x] && IGtQ[(n - 1)/2, 0]

Rule 3598

Int[sin[(e_.) + (f_.)*(x_)]^(m_.)*((a_) + (b_.)*tan[(e_.) + (f_.)*(x_)])^(n_.), x_Symbol] :> Int[Expand[Sin[e
+ f*x]^m*(a + b*Tan[e + f*x])^n, x], x] /; FreeQ[{a, b, e, f}, x] && IntegerQ[(m - 1)/2] && IGtQ[n, 0]

Rubi steps \begin{align*} \text {integral}& = \int \left (a \sin ^5(c+d x)+b \sin ^5(c+d x) \tan (c+d x)\right ) \, dx \\ & = a \int \sin ^5(c+d x) \, dx+b \int \sin ^5(c+d x) \tan (c+d x) \, dx \\ & = -\frac {a \text {Subst}\left (\int \left (1-2 x^2+x^4\right ) \, dx,x,\cos (c+d x)\right )}{d}+\frac {b \text {Subst}\left (\int \frac {x^6}{1-x^2} \, dx,x,\sin (c+d x)\right )}{d} \\ & = -\frac {a \cos (c+d x)}{d}+\frac {2 a \cos ^3(c+d x)}{3 d}-\frac {a \cos ^5(c+d x)}{5 d}+\frac {b \text {Subst}\left (\int \left (-1-x^2-x^4+\frac {1}{1-x^2}\right ) \, dx,x,\sin (c+d x)\right )}{d} \\ & = -\frac {a \cos (c+d x)}{d}+\frac {2 a \cos ^3(c+d x)}{3 d}-\frac {a \cos ^5(c+d x)}{5 d}-\frac {b \sin (c+d x)}{d}-\frac {b \sin ^3(c+d x)}{3 d}-\frac {b \sin ^5(c+d x)}{5 d}+\frac {b \text {Subst}\left (\int \frac {1}{1-x^2} \, dx,x,\sin (c+d x)\right )}{d} \\ & = \frac {b \text {arctanh}(\sin (c+d x))}{d}-\frac {a \cos (c+d x)}{d}+\frac {2 a \cos ^3(c+d x)}{3 d}-\frac {a \cos ^5(c+d x)}{5 d}-\frac {b \sin (c+d x)}{d}-\frac {b \sin ^3(c+d x)}{3 d}-\frac {b \sin ^5(c+d x)}{5 d} \\ \end{align*}

Mathematica [A] (verified)

Time = 0.09 (sec) , antiderivative size = 103, normalized size of antiderivative = 1.02 \[ \int \sin ^5(c+d x) (a+b \tan (c+d x)) \, dx=\frac {b \text {arctanh}(\sin (c+d x))}{d}-\frac {5 a \cos (c+d x)}{8 d}+\frac {5 a \cos (3 (c+d x))}{48 d}-\frac {a \cos (5 (c+d x))}{80 d}-\frac {b \sin (c+d x)}{d}-\frac {b \sin ^3(c+d x)}{3 d}-\frac {b \sin ^5(c+d x)}{5 d} \]

[In]

Integrate[Sin[c + d*x]^5*(a + b*Tan[c + d*x]),x]

[Out]

(b*ArcTanh[Sin[c + d*x]])/d - (5*a*Cos[c + d*x])/(8*d) + (5*a*Cos[3*(c + d*x)])/(48*d) - (a*Cos[5*(c + d*x)])/
(80*d) - (b*Sin[c + d*x])/d - (b*Sin[c + d*x]^3)/(3*d) - (b*Sin[c + d*x]^5)/(5*d)

Maple [A] (verified)

Time = 3.13 (sec) , antiderivative size = 80, normalized size of antiderivative = 0.79

method result size
derivativedivides \(\frac {-\frac {a \left (\frac {8}{3}+\sin ^{4}\left (d x +c \right )+\frac {4 \left (\sin ^{2}\left (d x +c \right )\right )}{3}\right ) \cos \left (d x +c \right )}{5}+b \left (-\frac {\left (\sin ^{5}\left (d x +c \right )\right )}{5}-\frac {\left (\sin ^{3}\left (d x +c \right )\right )}{3}-\sin \left (d x +c \right )+\ln \left (\sec \left (d x +c \right )+\tan \left (d x +c \right )\right )\right )}{d}\) \(80\)
default \(\frac {-\frac {a \left (\frac {8}{3}+\sin ^{4}\left (d x +c \right )+\frac {4 \left (\sin ^{2}\left (d x +c \right )\right )}{3}\right ) \cos \left (d x +c \right )}{5}+b \left (-\frac {\left (\sin ^{5}\left (d x +c \right )\right )}{5}-\frac {\left (\sin ^{3}\left (d x +c \right )\right )}{3}-\sin \left (d x +c \right )+\ln \left (\sec \left (d x +c \right )+\tan \left (d x +c \right )\right )\right )}{d}\) \(80\)
risch \(\frac {11 i {\mathrm e}^{i \left (d x +c \right )} b}{16 d}-\frac {5 \,{\mathrm e}^{i \left (d x +c \right )} a}{16 d}-\frac {11 i {\mathrm e}^{-i \left (d x +c \right )} b}{16 d}-\frac {5 \,{\mathrm e}^{-i \left (d x +c \right )} a}{16 d}+\frac {b \ln \left ({\mathrm e}^{i \left (d x +c \right )}+i\right )}{d}-\frac {b \ln \left ({\mathrm e}^{i \left (d x +c \right )}-i\right )}{d}-\frac {a \cos \left (5 d x +5 c \right )}{80 d}-\frac {b \sin \left (5 d x +5 c \right )}{80 d}+\frac {5 a \cos \left (3 d x +3 c \right )}{48 d}+\frac {7 b \sin \left (3 d x +3 c \right )}{48 d}\) \(161\)

[In]

int(sin(d*x+c)^5*(a+b*tan(d*x+c)),x,method=_RETURNVERBOSE)

[Out]

1/d*(-1/5*a*(8/3+sin(d*x+c)^4+4/3*sin(d*x+c)^2)*cos(d*x+c)+b*(-1/5*sin(d*x+c)^5-1/3*sin(d*x+c)^3-sin(d*x+c)+ln
(sec(d*x+c)+tan(d*x+c))))

Fricas [A] (verification not implemented)

none

Time = 0.27 (sec) , antiderivative size = 97, normalized size of antiderivative = 0.96 \[ \int \sin ^5(c+d x) (a+b \tan (c+d x)) \, dx=-\frac {6 \, a \cos \left (d x + c\right )^{5} - 20 \, a \cos \left (d x + c\right )^{3} + 30 \, a \cos \left (d x + c\right ) - 15 \, b \log \left (\sin \left (d x + c\right ) + 1\right ) + 15 \, b \log \left (-\sin \left (d x + c\right ) + 1\right ) + 2 \, {\left (3 \, b \cos \left (d x + c\right )^{4} - 11 \, b \cos \left (d x + c\right )^{2} + 23 \, b\right )} \sin \left (d x + c\right )}{30 \, d} \]

[In]

integrate(sin(d*x+c)^5*(a+b*tan(d*x+c)),x, algorithm="fricas")

[Out]

-1/30*(6*a*cos(d*x + c)^5 - 20*a*cos(d*x + c)^3 + 30*a*cos(d*x + c) - 15*b*log(sin(d*x + c) + 1) + 15*b*log(-s
in(d*x + c) + 1) + 2*(3*b*cos(d*x + c)^4 - 11*b*cos(d*x + c)^2 + 23*b)*sin(d*x + c))/d

Sympy [F]

\[ \int \sin ^5(c+d x) (a+b \tan (c+d x)) \, dx=\int \left (a + b \tan {\left (c + d x \right )}\right ) \sin ^{5}{\left (c + d x \right )}\, dx \]

[In]

integrate(sin(d*x+c)**5*(a+b*tan(d*x+c)),x)

[Out]

Integral((a + b*tan(c + d*x))*sin(c + d*x)**5, x)

Maxima [A] (verification not implemented)

none

Time = 0.21 (sec) , antiderivative size = 91, normalized size of antiderivative = 0.90 \[ \int \sin ^5(c+d x) (a+b \tan (c+d x)) \, dx=-\frac {2 \, {\left (3 \, \cos \left (d x + c\right )^{5} - 10 \, \cos \left (d x + c\right )^{3} + 15 \, \cos \left (d x + c\right )\right )} a + {\left (6 \, \sin \left (d x + c\right )^{5} + 10 \, \sin \left (d x + c\right )^{3} - 15 \, \log \left (\sin \left (d x + c\right ) + 1\right ) + 15 \, \log \left (\sin \left (d x + c\right ) - 1\right ) + 30 \, \sin \left (d x + c\right )\right )} b}{30 \, d} \]

[In]

integrate(sin(d*x+c)^5*(a+b*tan(d*x+c)),x, algorithm="maxima")

[Out]

-1/30*(2*(3*cos(d*x + c)^5 - 10*cos(d*x + c)^3 + 15*cos(d*x + c))*a + (6*sin(d*x + c)^5 + 10*sin(d*x + c)^3 -
15*log(sin(d*x + c) + 1) + 15*log(sin(d*x + c) - 1) + 30*sin(d*x + c))*b)/d

Giac [B] (verification not implemented)

Leaf count of result is larger than twice the leaf count of optimal. 10412 vs. \(2 (93) = 186\).

Time = 1.73 (sec) , antiderivative size = 10412, normalized size of antiderivative = 103.09 \[ \int \sin ^5(c+d x) (a+b \tan (c+d x)) \, dx=\text {Too large to display} \]

[In]

integrate(sin(d*x+c)^5*(a+b*tan(d*x+c)),x, algorithm="giac")

[Out]

-1/30*(15*b*log(2*(tan(1/2*d*x)^2*tan(1/2*c)^2 + 2*tan(1/2*d*x)^2*tan(1/2*c) + 2*tan(1/2*d*x)*tan(1/2*c)^2 + t
an(1/2*d*x)^2 + tan(1/2*c)^2 - 2*tan(1/2*d*x) - 2*tan(1/2*c) + 1)/(tan(1/2*d*x)^2*tan(1/2*c)^2 + tan(1/2*d*x)^
2 + tan(1/2*c)^2 + 1))*tan(1/2*d*x)^10*tan(1/2*c)^10 - 15*b*log(2*(tan(1/2*d*x)^2*tan(1/2*c)^2 - 2*tan(1/2*d*x
)^2*tan(1/2*c) - 2*tan(1/2*d*x)*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 2*tan(1/2*d*x) + 2*tan(1/2*c) +
 1)/(tan(1/2*d*x)^2*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 1))*tan(1/2*d*x)^10*tan(1/2*c)^10 + 16*a*ta
n(1/2*d*x)^10*tan(1/2*c)^10 + 75*b*log(2*(tan(1/2*d*x)^2*tan(1/2*c)^2 + 2*tan(1/2*d*x)^2*tan(1/2*c) + 2*tan(1/
2*d*x)*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 - 2*tan(1/2*d*x) - 2*tan(1/2*c) + 1)/(tan(1/2*d*x)^2*tan(1
/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 1))*tan(1/2*d*x)^10*tan(1/2*c)^8 - 75*b*log(2*(tan(1/2*d*x)^2*tan(1/
2*c)^2 - 2*tan(1/2*d*x)^2*tan(1/2*c) - 2*tan(1/2*d*x)*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 2*tan(1/2
*d*x) + 2*tan(1/2*c) + 1)/(tan(1/2*d*x)^2*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 1))*tan(1/2*d*x)^10*t
an(1/2*c)^8 - 60*b*tan(1/2*d*x)^10*tan(1/2*c)^9 + 75*b*log(2*(tan(1/2*d*x)^2*tan(1/2*c)^2 + 2*tan(1/2*d*x)^2*t
an(1/2*c) + 2*tan(1/2*d*x)*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 - 2*tan(1/2*d*x) - 2*tan(1/2*c) + 1)/(
tan(1/2*d*x)^2*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 1))*tan(1/2*d*x)^8*tan(1/2*c)^10 - 75*b*log(2*(t
an(1/2*d*x)^2*tan(1/2*c)^2 - 2*tan(1/2*d*x)^2*tan(1/2*c) - 2*tan(1/2*d*x)*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(
1/2*c)^2 + 2*tan(1/2*d*x) + 2*tan(1/2*c) + 1)/(tan(1/2*d*x)^2*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 1
))*tan(1/2*d*x)^8*tan(1/2*c)^10 - 60*b*tan(1/2*d*x)^9*tan(1/2*c)^10 + 80*a*tan(1/2*d*x)^10*tan(1/2*c)^8 + 80*a
*tan(1/2*d*x)^8*tan(1/2*c)^10 + 150*b*log(2*(tan(1/2*d*x)^2*tan(1/2*c)^2 + 2*tan(1/2*d*x)^2*tan(1/2*c) + 2*tan
(1/2*d*x)*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 - 2*tan(1/2*d*x) - 2*tan(1/2*c) + 1)/(tan(1/2*d*x)^2*ta
n(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 1))*tan(1/2*d*x)^10*tan(1/2*c)^6 - 150*b*log(2*(tan(1/2*d*x)^2*ta
n(1/2*c)^2 - 2*tan(1/2*d*x)^2*tan(1/2*c) - 2*tan(1/2*d*x)*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 2*tan
(1/2*d*x) + 2*tan(1/2*c) + 1)/(tan(1/2*d*x)^2*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 1))*tan(1/2*d*x)^
10*tan(1/2*c)^6 - 320*b*tan(1/2*d*x)^10*tan(1/2*c)^7 + 375*b*log(2*(tan(1/2*d*x)^2*tan(1/2*c)^2 + 2*tan(1/2*d*
x)^2*tan(1/2*c) + 2*tan(1/2*d*x)*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 - 2*tan(1/2*d*x) - 2*tan(1/2*c)
+ 1)/(tan(1/2*d*x)^2*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 1))*tan(1/2*d*x)^8*tan(1/2*c)^8 - 375*b*lo
g(2*(tan(1/2*d*x)^2*tan(1/2*c)^2 - 2*tan(1/2*d*x)^2*tan(1/2*c) - 2*tan(1/2*d*x)*tan(1/2*c)^2 + tan(1/2*d*x)^2
+ tan(1/2*c)^2 + 2*tan(1/2*d*x) + 2*tan(1/2*c) + 1)/(tan(1/2*d*x)^2*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)
^2 + 1))*tan(1/2*d*x)^8*tan(1/2*c)^8 - 420*b*tan(1/2*d*x)^9*tan(1/2*c)^8 - 420*b*tan(1/2*d*x)^8*tan(1/2*c)^9 +
 150*b*log(2*(tan(1/2*d*x)^2*tan(1/2*c)^2 + 2*tan(1/2*d*x)^2*tan(1/2*c) + 2*tan(1/2*d*x)*tan(1/2*c)^2 + tan(1/
2*d*x)^2 + tan(1/2*c)^2 - 2*tan(1/2*d*x) - 2*tan(1/2*c) + 1)/(tan(1/2*d*x)^2*tan(1/2*c)^2 + tan(1/2*d*x)^2 + t
an(1/2*c)^2 + 1))*tan(1/2*d*x)^6*tan(1/2*c)^10 - 150*b*log(2*(tan(1/2*d*x)^2*tan(1/2*c)^2 - 2*tan(1/2*d*x)^2*t
an(1/2*c) - 2*tan(1/2*d*x)*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 2*tan(1/2*d*x) + 2*tan(1/2*c) + 1)/(
tan(1/2*d*x)^2*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 1))*tan(1/2*d*x)^6*tan(1/2*c)^10 - 320*b*tan(1/2
*d*x)^7*tan(1/2*c)^10 + 160*a*tan(1/2*d*x)^10*tan(1/2*c)^6 + 400*a*tan(1/2*d*x)^8*tan(1/2*c)^8 + 160*a*tan(1/2
*d*x)^6*tan(1/2*c)^10 + 150*b*log(2*(tan(1/2*d*x)^2*tan(1/2*c)^2 + 2*tan(1/2*d*x)^2*tan(1/2*c) + 2*tan(1/2*d*x
)*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 - 2*tan(1/2*d*x) - 2*tan(1/2*c) + 1)/(tan(1/2*d*x)^2*tan(1/2*c)
^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 1))*tan(1/2*d*x)^10*tan(1/2*c)^4 - 150*b*log(2*(tan(1/2*d*x)^2*tan(1/2*c)
^2 - 2*tan(1/2*d*x)^2*tan(1/2*c) - 2*tan(1/2*d*x)*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 2*tan(1/2*d*x
) + 2*tan(1/2*c) + 1)/(tan(1/2*d*x)^2*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 1))*tan(1/2*d*x)^10*tan(1
/2*c)^4 - 712*b*tan(1/2*d*x)^10*tan(1/2*c)^5 + 750*b*log(2*(tan(1/2*d*x)^2*tan(1/2*c)^2 + 2*tan(1/2*d*x)^2*tan
(1/2*c) + 2*tan(1/2*d*x)*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 - 2*tan(1/2*d*x) - 2*tan(1/2*c) + 1)/(ta
n(1/2*d*x)^2*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 1))*tan(1/2*d*x)^8*tan(1/2*c)^6 - 750*b*log(2*(tan
(1/2*d*x)^2*tan(1/2*c)^2 - 2*tan(1/2*d*x)^2*tan(1/2*c) - 2*tan(1/2*d*x)*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/
2*c)^2 + 2*tan(1/2*d*x) + 2*tan(1/2*c) + 1)/(tan(1/2*d*x)^2*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 1))
*tan(1/2*d*x)^8*tan(1/2*c)^6 - 1320*b*tan(1/2*d*x)^9*tan(1/2*c)^6 - 2560*b*tan(1/2*d*x)^8*tan(1/2*c)^7 + 750*b
*log(2*(tan(1/2*d*x)^2*tan(1/2*c)^2 + 2*tan(1/2*d*x)^2*tan(1/2*c) + 2*tan(1/2*d*x)*tan(1/2*c)^2 + tan(1/2*d*x)
^2 + tan(1/2*c)^2 - 2*tan(1/2*d*x) - 2*tan(1/2*c) + 1)/(tan(1/2*d*x)^2*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2
*c)^2 + 1))*tan(1/2*d*x)^6*tan(1/2*c)^8 - 750*b*log(2*(tan(1/2*d*x)^2*tan(1/2*c)^2 - 2*tan(1/2*d*x)^2*tan(1/2*
c) - 2*tan(1/2*d*x)*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 2*tan(1/2*d*x) + 2*tan(1/2*c) + 1)/(tan(1/2
*d*x)^2*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 1))*tan(1/2*d*x)^6*tan(1/2*c)^8 - 2560*b*tan(1/2*d*x)^7
*tan(1/2*c)^8 - 1320*b*tan(1/2*d*x)^6*tan(1/2*c)^9 + 150*b*log(2*(tan(1/2*d*x)^2*tan(1/2*c)^2 + 2*tan(1/2*d*x)
^2*tan(1/2*c) + 2*tan(1/2*d*x)*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 - 2*tan(1/2*d*x) - 2*tan(1/2*c) +
1)/(tan(1/2*d*x)^2*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 1))*tan(1/2*d*x)^4*tan(1/2*c)^10 - 150*b*log
(2*(tan(1/2*d*x)^2*tan(1/2*c)^2 - 2*tan(1/2*d*x)^2*tan(1/2*c) - 2*tan(1/2*d*x)*tan(1/2*c)^2 + tan(1/2*d*x)^2 +
 tan(1/2*c)^2 + 2*tan(1/2*d*x) + 2*tan(1/2*c) + 1)/(tan(1/2*d*x)^2*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^
2 + 1))*tan(1/2*d*x)^4*tan(1/2*c)^10 - 712*b*tan(1/2*d*x)^5*tan(1/2*c)^10 - 160*a*tan(1/2*d*x)^10*tan(1/2*c)^4
 - 1920*a*tan(1/2*d*x)^9*tan(1/2*c)^5 - 4000*a*tan(1/2*d*x)^8*tan(1/2*c)^6 - 6400*a*tan(1/2*d*x)^7*tan(1/2*c)^
7 - 4000*a*tan(1/2*d*x)^6*tan(1/2*c)^8 - 1920*a*tan(1/2*d*x)^5*tan(1/2*c)^9 - 160*a*tan(1/2*d*x)^4*tan(1/2*c)^
10 + 75*b*log(2*(tan(1/2*d*x)^2*tan(1/2*c)^2 + 2*tan(1/2*d*x)^2*tan(1/2*c) + 2*tan(1/2*d*x)*tan(1/2*c)^2 + tan
(1/2*d*x)^2 + tan(1/2*c)^2 - 2*tan(1/2*d*x) - 2*tan(1/2*c) + 1)/(tan(1/2*d*x)^2*tan(1/2*c)^2 + tan(1/2*d*x)^2
+ tan(1/2*c)^2 + 1))*tan(1/2*d*x)^10*tan(1/2*c)^2 - 75*b*log(2*(tan(1/2*d*x)^2*tan(1/2*c)^2 - 2*tan(1/2*d*x)^2
*tan(1/2*c) - 2*tan(1/2*d*x)*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 2*tan(1/2*d*x) + 2*tan(1/2*c) + 1)
/(tan(1/2*d*x)^2*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 1))*tan(1/2*d*x)^10*tan(1/2*c)^2 - 320*b*tan(1
/2*d*x)^10*tan(1/2*c)^3 + 750*b*log(2*(tan(1/2*d*x)^2*tan(1/2*c)^2 + 2*tan(1/2*d*x)^2*tan(1/2*c) + 2*tan(1/2*d
*x)*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 - 2*tan(1/2*d*x) - 2*tan(1/2*c) + 1)/(tan(1/2*d*x)^2*tan(1/2*
c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 1))*tan(1/2*d*x)^8*tan(1/2*c)^4 - 750*b*log(2*(tan(1/2*d*x)^2*tan(1/2*c
)^2 - 2*tan(1/2*d*x)^2*tan(1/2*c) - 2*tan(1/2*d*x)*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 2*tan(1/2*d*
x) + 2*tan(1/2*c) + 1)/(tan(1/2*d*x)^2*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 1))*tan(1/2*d*x)^8*tan(1
/2*c)^4 + 1320*b*tan(1/2*d*x)^9*tan(1/2*c)^4 + 4360*b*tan(1/2*d*x)^8*tan(1/2*c)^5 + 1500*b*log(2*(tan(1/2*d*x)
^2*tan(1/2*c)^2 + 2*tan(1/2*d*x)^2*tan(1/2*c) + 2*tan(1/2*d*x)*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 -
2*tan(1/2*d*x) - 2*tan(1/2*c) + 1)/(tan(1/2*d*x)^2*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 1))*tan(1/2*
d*x)^6*tan(1/2*c)^6 - 1500*b*log(2*(tan(1/2*d*x)^2*tan(1/2*c)^2 - 2*tan(1/2*d*x)^2*tan(1/2*c) - 2*tan(1/2*d*x)
*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 2*tan(1/2*d*x) + 2*tan(1/2*c) + 1)/(tan(1/2*d*x)^2*tan(1/2*c)^
2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 1))*tan(1/2*d*x)^6*tan(1/2*c)^6 + 9280*b*tan(1/2*d*x)^7*tan(1/2*c)^6 + 928
0*b*tan(1/2*d*x)^6*tan(1/2*c)^7 + 750*b*log(2*(tan(1/2*d*x)^2*tan(1/2*c)^2 + 2*tan(1/2*d*x)^2*tan(1/2*c) + 2*t
an(1/2*d*x)*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 - 2*tan(1/2*d*x) - 2*tan(1/2*c) + 1)/(tan(1/2*d*x)^2*
tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 1))*tan(1/2*d*x)^4*tan(1/2*c)^8 - 750*b*log(2*(tan(1/2*d*x)^2*t
an(1/2*c)^2 - 2*tan(1/2*d*x)^2*tan(1/2*c) - 2*tan(1/2*d*x)*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 2*ta
n(1/2*d*x) + 2*tan(1/2*c) + 1)/(tan(1/2*d*x)^2*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 1))*tan(1/2*d*x)
^4*tan(1/2*c)^8 + 4360*b*tan(1/2*d*x)^5*tan(1/2*c)^8 + 1320*b*tan(1/2*d*x)^4*tan(1/2*c)^9 + 75*b*log(2*(tan(1/
2*d*x)^2*tan(1/2*c)^2 + 2*tan(1/2*d*x)^2*tan(1/2*c) + 2*tan(1/2*d*x)*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c
)^2 - 2*tan(1/2*d*x) - 2*tan(1/2*c) + 1)/(tan(1/2*d*x)^2*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 1))*ta
n(1/2*d*x)^2*tan(1/2*c)^10 - 75*b*log(2*(tan(1/2*d*x)^2*tan(1/2*c)^2 - 2*tan(1/2*d*x)^2*tan(1/2*c) - 2*tan(1/2
*d*x)*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 2*tan(1/2*d*x) + 2*tan(1/2*c) + 1)/(tan(1/2*d*x)^2*tan(1/
2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 1))*tan(1/2*d*x)^2*tan(1/2*c)^10 - 320*b*tan(1/2*d*x)^3*tan(1/2*c)^10
 - 80*a*tan(1/2*d*x)^10*tan(1/2*c)^2 + 4000*a*tan(1/2*d*x)^8*tan(1/2*c)^4 + 10240*a*tan(1/2*d*x)^7*tan(1/2*c)^
5 + 16000*a*tan(1/2*d*x)^6*tan(1/2*c)^6 + 10240*a*tan(1/2*d*x)^5*tan(1/2*c)^7 + 4000*a*tan(1/2*d*x)^4*tan(1/2*
c)^8 - 80*a*tan(1/2*d*x)^2*tan(1/2*c)^10 + 15*b*log(2*(tan(1/2*d*x)^2*tan(1/2*c)^2 + 2*tan(1/2*d*x)^2*tan(1/2*
c) + 2*tan(1/2*d*x)*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 - 2*tan(1/2*d*x) - 2*tan(1/2*c) + 1)/(tan(1/2
*d*x)^2*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 1))*tan(1/2*d*x)^10 - 15*b*log(2*(tan(1/2*d*x)^2*tan(1/
2*c)^2 - 2*tan(1/2*d*x)^2*tan(1/2*c) - 2*tan(1/2*d*x)*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 2*tan(1/2
*d*x) + 2*tan(1/2*c) + 1)/(tan(1/2*d*x)^2*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 1))*tan(1/2*d*x)^10 -
 60*b*tan(1/2*d*x)^10*tan(1/2*c) + 375*b*log(2*(tan(1/2*d*x)^2*tan(1/2*c)^2 + 2*tan(1/2*d*x)^2*tan(1/2*c) + 2*
tan(1/2*d*x)*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 - 2*tan(1/2*d*x) - 2*tan(1/2*c) + 1)/(tan(1/2*d*x)^2
*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 1))*tan(1/2*d*x)^8*tan(1/2*c)^2 - 375*b*log(2*(tan(1/2*d*x)^2*
tan(1/2*c)^2 - 2*tan(1/2*d*x)^2*tan(1/2*c) - 2*tan(1/2*d*x)*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 2*t
an(1/2*d*x) + 2*tan(1/2*c) + 1)/(tan(1/2*d*x)^2*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 1))*tan(1/2*d*x
)^8*tan(1/2*c)^2 + 420*b*tan(1/2*d*x)^9*tan(1/2*c)^2 - 2560*b*tan(1/2*d*x)^8*tan(1/2*c)^3 + 1500*b*log(2*(tan(
1/2*d*x)^2*tan(1/2*c)^2 + 2*tan(1/2*d*x)^2*tan(1/2*c) + 2*tan(1/2*d*x)*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2
*c)^2 - 2*tan(1/2*d*x) - 2*tan(1/2*c) + 1)/(tan(1/2*d*x)^2*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 1))*
tan(1/2*d*x)^6*tan(1/2*c)^4 - 1500*b*log(2*(tan(1/2*d*x)^2*tan(1/2*c)^2 - 2*tan(1/2*d*x)^2*tan(1/2*c) - 2*tan(
1/2*d*x)*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 2*tan(1/2*d*x) + 2*tan(1/2*c) + 1)/(tan(1/2*d*x)^2*tan
(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 1))*tan(1/2*d*x)^6*tan(1/2*c)^4 - 9280*b*tan(1/2*d*x)^7*tan(1/2*c)
^4 - 19120*b*tan(1/2*d*x)^6*tan(1/2*c)^5 + 1500*b*log(2*(tan(1/2*d*x)^2*tan(1/2*c)^2 + 2*tan(1/2*d*x)^2*tan(1/
2*c) + 2*tan(1/2*d*x)*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 - 2*tan(1/2*d*x) - 2*tan(1/2*c) + 1)/(tan(1
/2*d*x)^2*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 1))*tan(1/2*d*x)^4*tan(1/2*c)^6 - 1500*b*log(2*(tan(1
/2*d*x)^2*tan(1/2*c)^2 - 2*tan(1/2*d*x)^2*tan(1/2*c) - 2*tan(1/2*d*x)*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*
c)^2 + 2*tan(1/2*d*x) + 2*tan(1/2*c) + 1)/(tan(1/2*d*x)^2*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 1))*t
an(1/2*d*x)^4*tan(1/2*c)^6 - 19120*b*tan(1/2*d*x)^5*tan(1/2*c)^6 - 9280*b*tan(1/2*d*x)^4*tan(1/2*c)^7 + 375*b*
log(2*(tan(1/2*d*x)^2*tan(1/2*c)^2 + 2*tan(1/2*d*x)^2*tan(1/2*c) + 2*tan(1/2*d*x)*tan(1/2*c)^2 + tan(1/2*d*x)^
2 + tan(1/2*c)^2 - 2*tan(1/2*d*x) - 2*tan(1/2*c) + 1)/(tan(1/2*d*x)^2*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*
c)^2 + 1))*tan(1/2*d*x)^2*tan(1/2*c)^8 - 375*b*log(2*(tan(1/2*d*x)^2*tan(1/2*c)^2 - 2*tan(1/2*d*x)^2*tan(1/2*c
) - 2*tan(1/2*d*x)*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 2*tan(1/2*d*x) + 2*tan(1/2*c) + 1)/(tan(1/2*
d*x)^2*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 1))*tan(1/2*d*x)^2*tan(1/2*c)^8 - 2560*b*tan(1/2*d*x)^3*
tan(1/2*c)^8 + 420*b*tan(1/2*d*x)^2*tan(1/2*c)^9 + 15*b*log(2*(tan(1/2*d*x)^2*tan(1/2*c)^2 + 2*tan(1/2*d*x)^2*
tan(1/2*c) + 2*tan(1/2*d*x)*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 - 2*tan(1/2*d*x) - 2*tan(1/2*c) + 1)/
(tan(1/2*d*x)^2*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 1))*tan(1/2*c)^10 - 15*b*log(2*(tan(1/2*d*x)^2*
tan(1/2*c)^2 - 2*tan(1/2*d*x)^2*tan(1/2*c) - 2*tan(1/2*d*x)*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 2*t
an(1/2*d*x) + 2*tan(1/2*c) + 1)/(tan(1/2*d*x)^2*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 1))*tan(1/2*c)^
10 - 60*b*tan(1/2*d*x)*tan(1/2*c)^10 - 16*a*tan(1/2*d*x)^10 - 400*a*tan(1/2*d*x)^8*tan(1/2*c)^2 - 6400*a*tan(1
/2*d*x)^7*tan(1/2*c)^3 - 16000*a*tan(1/2*d*x)^6*tan(1/2*c)^4 - 24064*a*tan(1/2*d*x)^5*tan(1/2*c)^5 - 16000*a*t
an(1/2*d*x)^4*tan(1/2*c)^6 - 6400*a*tan(1/2*d*x)^3*tan(1/2*c)^7 - 400*a*tan(1/2*d*x)^2*tan(1/2*c)^8 - 16*a*tan
(1/2*c)^10 + 75*b*log(2*(tan(1/2*d*x)^2*tan(1/2*c)^2 + 2*tan(1/2*d*x)^2*tan(1/2*c) + 2*tan(1/2*d*x)*tan(1/2*c)
^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 - 2*tan(1/2*d*x) - 2*tan(1/2*c) + 1)/(tan(1/2*d*x)^2*tan(1/2*c)^2 + tan(1/2
*d*x)^2 + tan(1/2*c)^2 + 1))*tan(1/2*d*x)^8 - 75*b*log(2*(tan(1/2*d*x)^2*tan(1/2*c)^2 - 2*tan(1/2*d*x)^2*tan(1
/2*c) - 2*tan(1/2*d*x)*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 2*tan(1/2*d*x) + 2*tan(1/2*c) + 1)/(tan(
1/2*d*x)^2*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 1))*tan(1/2*d*x)^8 + 60*b*tan(1/2*d*x)^9 - 420*b*tan
(1/2*d*x)^8*tan(1/2*c) + 750*b*log(2*(tan(1/2*d*x)^2*tan(1/2*c)^2 + 2*tan(1/2*d*x)^2*tan(1/2*c) + 2*tan(1/2*d*
x)*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 - 2*tan(1/2*d*x) - 2*tan(1/2*c) + 1)/(tan(1/2*d*x)^2*tan(1/2*c
)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 1))*tan(1/2*d*x)^6*tan(1/2*c)^2 - 750*b*log(2*(tan(1/2*d*x)^2*tan(1/2*c)
^2 - 2*tan(1/2*d*x)^2*tan(1/2*c) - 2*tan(1/2*d*x)*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 2*tan(1/2*d*x
) + 2*tan(1/2*c) + 1)/(tan(1/2*d*x)^2*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 1))*tan(1/2*d*x)^6*tan(1/
2*c)^2 + 2560*b*tan(1/2*d*x)^7*tan(1/2*c)^2 + 9280*b*tan(1/2*d*x)^6*tan(1/2*c)^3 + 1500*b*log(2*(tan(1/2*d*x)^
2*tan(1/2*c)^2 + 2*tan(1/2*d*x)^2*tan(1/2*c) + 2*tan(1/2*d*x)*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 - 2
*tan(1/2*d*x) - 2*tan(1/2*c) + 1)/(tan(1/2*d*x)^2*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 1))*tan(1/2*d
*x)^4*tan(1/2*c)^4 - 1500*b*log(2*(tan(1/2*d*x)^2*tan(1/2*c)^2 - 2*tan(1/2*d*x)^2*tan(1/2*c) - 2*tan(1/2*d*x)*
tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 2*tan(1/2*d*x) + 2*tan(1/2*c) + 1)/(tan(1/2*d*x)^2*tan(1/2*c)^2
 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 1))*tan(1/2*d*x)^4*tan(1/2*c)^4 + 19120*b*tan(1/2*d*x)^5*tan(1/2*c)^4 + 191
20*b*tan(1/2*d*x)^4*tan(1/2*c)^5 + 750*b*log(2*(tan(1/2*d*x)^2*tan(1/2*c)^2 + 2*tan(1/2*d*x)^2*tan(1/2*c) + 2*
tan(1/2*d*x)*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 - 2*tan(1/2*d*x) - 2*tan(1/2*c) + 1)/(tan(1/2*d*x)^2
*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 1))*tan(1/2*d*x)^2*tan(1/2*c)^6 - 750*b*log(2*(tan(1/2*d*x)^2*
tan(1/2*c)^2 - 2*tan(1/2*d*x)^2*tan(1/2*c) - 2*tan(1/2*d*x)*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 2*t
an(1/2*d*x) + 2*tan(1/2*c) + 1)/(tan(1/2*d*x)^2*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 1))*tan(1/2*d*x
)^2*tan(1/2*c)^6 + 9280*b*tan(1/2*d*x)^3*tan(1/2*c)^6 + 2560*b*tan(1/2*d*x)^2*tan(1/2*c)^7 + 75*b*log(2*(tan(1
/2*d*x)^2*tan(1/2*c)^2 + 2*tan(1/2*d*x)^2*tan(1/2*c) + 2*tan(1/2*d*x)*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*
c)^2 - 2*tan(1/2*d*x) - 2*tan(1/2*c) + 1)/(tan(1/2*d*x)^2*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 1))*t
an(1/2*c)^8 - 75*b*log(2*(tan(1/2*d*x)^2*tan(1/2*c)^2 - 2*tan(1/2*d*x)^2*tan(1/2*c) - 2*tan(1/2*d*x)*tan(1/2*c
)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 2*tan(1/2*d*x) + 2*tan(1/2*c) + 1)/(tan(1/2*d*x)^2*tan(1/2*c)^2 + tan(1/
2*d*x)^2 + tan(1/2*c)^2 + 1))*tan(1/2*c)^8 - 420*b*tan(1/2*d*x)*tan(1/2*c)^8 + 60*b*tan(1/2*c)^9 - 80*a*tan(1/
2*d*x)^8 + 4000*a*tan(1/2*d*x)^6*tan(1/2*c)^2 + 10240*a*tan(1/2*d*x)^5*tan(1/2*c)^3 + 16000*a*tan(1/2*d*x)^4*t
an(1/2*c)^4 + 10240*a*tan(1/2*d*x)^3*tan(1/2*c)^5 + 4000*a*tan(1/2*d*x)^2*tan(1/2*c)^6 - 80*a*tan(1/2*c)^8 + 1
50*b*log(2*(tan(1/2*d*x)^2*tan(1/2*c)^2 + 2*tan(1/2*d*x)^2*tan(1/2*c) + 2*tan(1/2*d*x)*tan(1/2*c)^2 + tan(1/2*
d*x)^2 + tan(1/2*c)^2 - 2*tan(1/2*d*x) - 2*tan(1/2*c) + 1)/(tan(1/2*d*x)^2*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan
(1/2*c)^2 + 1))*tan(1/2*d*x)^6 - 150*b*log(2*(tan(1/2*d*x)^2*tan(1/2*c)^2 - 2*tan(1/2*d*x)^2*tan(1/2*c) - 2*ta
n(1/2*d*x)*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 2*tan(1/2*d*x) + 2*tan(1/2*c) + 1)/(tan(1/2*d*x)^2*t
an(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 1))*tan(1/2*d*x)^6 + 320*b*tan(1/2*d*x)^7 - 1320*b*tan(1/2*d*x)^
6*tan(1/2*c) + 750*b*log(2*(tan(1/2*d*x)^2*tan(1/2*c)^2 + 2*tan(1/2*d*x)^2*tan(1/2*c) + 2*tan(1/2*d*x)*tan(1/2
*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 - 2*tan(1/2*d*x) - 2*tan(1/2*c) + 1)/(tan(1/2*d*x)^2*tan(1/2*c)^2 + tan(
1/2*d*x)^2 + tan(1/2*c)^2 + 1))*tan(1/2*d*x)^4*tan(1/2*c)^2 - 750*b*log(2*(tan(1/2*d*x)^2*tan(1/2*c)^2 - 2*tan
(1/2*d*x)^2*tan(1/2*c) - 2*tan(1/2*d*x)*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 2*tan(1/2*d*x) + 2*tan(
1/2*c) + 1)/(tan(1/2*d*x)^2*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 1))*tan(1/2*d*x)^4*tan(1/2*c)^2 - 4
360*b*tan(1/2*d*x)^5*tan(1/2*c)^2 - 9280*b*tan(1/2*d*x)^4*tan(1/2*c)^3 + 750*b*log(2*(tan(1/2*d*x)^2*tan(1/2*c
)^2 + 2*tan(1/2*d*x)^2*tan(1/2*c) + 2*tan(1/2*d*x)*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 - 2*tan(1/2*d*
x) - 2*tan(1/2*c) + 1)/(tan(1/2*d*x)^2*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 1))*tan(1/2*d*x)^2*tan(1
/2*c)^4 - 750*b*log(2*(tan(1/2*d*x)^2*tan(1/2*c)^2 - 2*tan(1/2*d*x)^2*tan(1/2*c) - 2*tan(1/2*d*x)*tan(1/2*c)^2
 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 2*tan(1/2*d*x) + 2*tan(1/2*c) + 1)/(tan(1/2*d*x)^2*tan(1/2*c)^2 + tan(1/2*d
*x)^2 + tan(1/2*c)^2 + 1))*tan(1/2*d*x)^2*tan(1/2*c)^4 - 9280*b*tan(1/2*d*x)^3*tan(1/2*c)^4 - 4360*b*tan(1/2*d
*x)^2*tan(1/2*c)^5 + 150*b*log(2*(tan(1/2*d*x)^2*tan(1/2*c)^2 + 2*tan(1/2*d*x)^2*tan(1/2*c) + 2*tan(1/2*d*x)*t
an(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 - 2*tan(1/2*d*x) - 2*tan(1/2*c) + 1)/(tan(1/2*d*x)^2*tan(1/2*c)^2
+ tan(1/2*d*x)^2 + tan(1/2*c)^2 + 1))*tan(1/2*c)^6 - 150*b*log(2*(tan(1/2*d*x)^2*tan(1/2*c)^2 - 2*tan(1/2*d*x)
^2*tan(1/2*c) - 2*tan(1/2*d*x)*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 2*tan(1/2*d*x) + 2*tan(1/2*c) +
1)/(tan(1/2*d*x)^2*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 1))*tan(1/2*c)^6 - 1320*b*tan(1/2*d*x)*tan(1
/2*c)^6 + 320*b*tan(1/2*c)^7 - 160*a*tan(1/2*d*x)^6 - 1920*a*tan(1/2*d*x)^5*tan(1/2*c) - 4000*a*tan(1/2*d*x)^4
*tan(1/2*c)^2 - 6400*a*tan(1/2*d*x)^3*tan(1/2*c)^3 - 4000*a*tan(1/2*d*x)^2*tan(1/2*c)^4 - 1920*a*tan(1/2*d*x)*
tan(1/2*c)^5 - 160*a*tan(1/2*c)^6 + 150*b*log(2*(tan(1/2*d*x)^2*tan(1/2*c)^2 + 2*tan(1/2*d*x)^2*tan(1/2*c) + 2
*tan(1/2*d*x)*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 - 2*tan(1/2*d*x) - 2*tan(1/2*c) + 1)/(tan(1/2*d*x)^
2*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 1))*tan(1/2*d*x)^4 - 150*b*log(2*(tan(1/2*d*x)^2*tan(1/2*c)^2
 - 2*tan(1/2*d*x)^2*tan(1/2*c) - 2*tan(1/2*d*x)*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 2*tan(1/2*d*x)
+ 2*tan(1/2*c) + 1)/(tan(1/2*d*x)^2*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 1))*tan(1/2*d*x)^4 + 712*b*
tan(1/2*d*x)^5 + 1320*b*tan(1/2*d*x)^4*tan(1/2*c) + 375*b*log(2*(tan(1/2*d*x)^2*tan(1/2*c)^2 + 2*tan(1/2*d*x)^
2*tan(1/2*c) + 2*tan(1/2*d*x)*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 - 2*tan(1/2*d*x) - 2*tan(1/2*c) + 1
)/(tan(1/2*d*x)^2*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 1))*tan(1/2*d*x)^2*tan(1/2*c)^2 - 375*b*log(2
*(tan(1/2*d*x)^2*tan(1/2*c)^2 - 2*tan(1/2*d*x)^2*tan(1/2*c) - 2*tan(1/2*d*x)*tan(1/2*c)^2 + tan(1/2*d*x)^2 + t
an(1/2*c)^2 + 2*tan(1/2*d*x) + 2*tan(1/2*c) + 1)/(tan(1/2*d*x)^2*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2
+ 1))*tan(1/2*d*x)^2*tan(1/2*c)^2 + 2560*b*tan(1/2*d*x)^3*tan(1/2*c)^2 + 2560*b*tan(1/2*d*x)^2*tan(1/2*c)^3 +
150*b*log(2*(tan(1/2*d*x)^2*tan(1/2*c)^2 + 2*tan(1/2*d*x)^2*tan(1/2*c) + 2*tan(1/2*d*x)*tan(1/2*c)^2 + tan(1/2
*d*x)^2 + tan(1/2*c)^2 - 2*tan(1/2*d*x) - 2*tan(1/2*c) + 1)/(tan(1/2*d*x)^2*tan(1/2*c)^2 + tan(1/2*d*x)^2 + ta
n(1/2*c)^2 + 1))*tan(1/2*c)^4 - 150*b*log(2*(tan(1/2*d*x)^2*tan(1/2*c)^2 - 2*tan(1/2*d*x)^2*tan(1/2*c) - 2*tan
(1/2*d*x)*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 2*tan(1/2*d*x) + 2*tan(1/2*c) + 1)/(tan(1/2*d*x)^2*ta
n(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 1))*tan(1/2*c)^4 + 1320*b*tan(1/2*d*x)*tan(1/2*c)^4 + 712*b*tan(1
/2*c)^5 + 160*a*tan(1/2*d*x)^4 + 400*a*tan(1/2*d*x)^2*tan(1/2*c)^2 + 160*a*tan(1/2*c)^4 + 75*b*log(2*(tan(1/2*
d*x)^2*tan(1/2*c)^2 + 2*tan(1/2*d*x)^2*tan(1/2*c) + 2*tan(1/2*d*x)*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^
2 - 2*tan(1/2*d*x) - 2*tan(1/2*c) + 1)/(tan(1/2*d*x)^2*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 1))*tan(
1/2*d*x)^2 - 75*b*log(2*(tan(1/2*d*x)^2*tan(1/2*c)^2 - 2*tan(1/2*d*x)^2*tan(1/2*c) - 2*tan(1/2*d*x)*tan(1/2*c)
^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 2*tan(1/2*d*x) + 2*tan(1/2*c) + 1)/(tan(1/2*d*x)^2*tan(1/2*c)^2 + tan(1/2
*d*x)^2 + tan(1/2*c)^2 + 1))*tan(1/2*d*x)^2 + 320*b*tan(1/2*d*x)^3 + 420*b*tan(1/2*d*x)^2*tan(1/2*c) + 75*b*lo
g(2*(tan(1/2*d*x)^2*tan(1/2*c)^2 + 2*tan(1/2*d*x)^2*tan(1/2*c) + 2*tan(1/2*d*x)*tan(1/2*c)^2 + tan(1/2*d*x)^2
+ tan(1/2*c)^2 - 2*tan(1/2*d*x) - 2*tan(1/2*c) + 1)/(tan(1/2*d*x)^2*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)
^2 + 1))*tan(1/2*c)^2 - 75*b*log(2*(tan(1/2*d*x)^2*tan(1/2*c)^2 - 2*tan(1/2*d*x)^2*tan(1/2*c) - 2*tan(1/2*d*x)
*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 2*tan(1/2*d*x) + 2*tan(1/2*c) + 1)/(tan(1/2*d*x)^2*tan(1/2*c)^
2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 1))*tan(1/2*c)^2 + 420*b*tan(1/2*d*x)*tan(1/2*c)^2 + 320*b*tan(1/2*c)^3 +
80*a*tan(1/2*d*x)^2 + 80*a*tan(1/2*c)^2 + 15*b*log(2*(tan(1/2*d*x)^2*tan(1/2*c)^2 + 2*tan(1/2*d*x)^2*tan(1/2*c
) + 2*tan(1/2*d*x)*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 - 2*tan(1/2*d*x) - 2*tan(1/2*c) + 1)/(tan(1/2*
d*x)^2*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 1)) - 15*b*log(2*(tan(1/2*d*x)^2*tan(1/2*c)^2 - 2*tan(1/
2*d*x)^2*tan(1/2*c) - 2*tan(1/2*d*x)*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 2*tan(1/2*d*x) + 2*tan(1/2
*c) + 1)/(tan(1/2*d*x)^2*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 1)) + 60*b*tan(1/2*d*x) + 60*b*tan(1/2
*c) + 16*a)/(d*tan(1/2*d*x)^10*tan(1/2*c)^10 + 5*d*tan(1/2*d*x)^10*tan(1/2*c)^8 + 5*d*tan(1/2*d*x)^8*tan(1/2*c
)^10 + 10*d*tan(1/2*d*x)^10*tan(1/2*c)^6 + 25*d*tan(1/2*d*x)^8*tan(1/2*c)^8 + 10*d*tan(1/2*d*x)^6*tan(1/2*c)^1
0 + 10*d*tan(1/2*d*x)^10*tan(1/2*c)^4 + 50*d*tan(1/2*d*x)^8*tan(1/2*c)^6 + 50*d*tan(1/2*d*x)^6*tan(1/2*c)^8 +
10*d*tan(1/2*d*x)^4*tan(1/2*c)^10 + 5*d*tan(1/2*d*x)^10*tan(1/2*c)^2 + 50*d*tan(1/2*d*x)^8*tan(1/2*c)^4 + 100*
d*tan(1/2*d*x)^6*tan(1/2*c)^6 + 50*d*tan(1/2*d*x)^4*tan(1/2*c)^8 + 5*d*tan(1/2*d*x)^2*tan(1/2*c)^10 + d*tan(1/
2*d*x)^10 + 25*d*tan(1/2*d*x)^8*tan(1/2*c)^2 + 100*d*tan(1/2*d*x)^6*tan(1/2*c)^4 + 100*d*tan(1/2*d*x)^4*tan(1/
2*c)^6 + 25*d*tan(1/2*d*x)^2*tan(1/2*c)^8 + d*tan(1/2*c)^10 + 5*d*tan(1/2*d*x)^8 + 50*d*tan(1/2*d*x)^6*tan(1/2
*c)^2 + 100*d*tan(1/2*d*x)^4*tan(1/2*c)^4 + 50*d*tan(1/2*d*x)^2*tan(1/2*c)^6 + 5*d*tan(1/2*c)^8 + 10*d*tan(1/2
*d*x)^6 + 50*d*tan(1/2*d*x)^4*tan(1/2*c)^2 + 50*d*tan(1/2*d*x)^2*tan(1/2*c)^4 + 10*d*tan(1/2*c)^6 + 10*d*tan(1
/2*d*x)^4 + 25*d*tan(1/2*d*x)^2*tan(1/2*c)^2 + 10*d*tan(1/2*c)^4 + 5*d*tan(1/2*d*x)^2 + 5*d*tan(1/2*c)^2 + d)

Mupad [B] (verification not implemented)

Time = 5.22 (sec) , antiderivative size = 121, normalized size of antiderivative = 1.20 \[ \int \sin ^5(c+d x) (a+b \tan (c+d x)) \, dx=\frac {2\,b\,\mathrm {atanh}\left (\frac {\sin \left (\frac {c}{2}+\frac {d\,x}{2}\right )}{\cos \left (\frac {c}{2}+\frac {d\,x}{2}\right )}\right )}{d}+\frac {2\,a\,{\cos \left (c+d\,x\right )}^3}{3\,d}-\frac {a\,{\cos \left (c+d\,x\right )}^5}{5\,d}-\frac {a\,\cos \left (c+d\,x\right )}{d}-\frac {23\,b\,\sin \left (c+d\,x\right )}{15\,d}+\frac {11\,b\,{\cos \left (c+d\,x\right )}^2\,\sin \left (c+d\,x\right )}{15\,d}-\frac {b\,{\cos \left (c+d\,x\right )}^4\,\sin \left (c+d\,x\right )}{5\,d} \]

[In]

int(sin(c + d*x)^5*(a + b*tan(c + d*x)),x)

[Out]

(2*b*atanh(sin(c/2 + (d*x)/2)/cos(c/2 + (d*x)/2)))/d + (2*a*cos(c + d*x)^3)/(3*d) - (a*cos(c + d*x)^5)/(5*d) -
 (a*cos(c + d*x))/d - (23*b*sin(c + d*x))/(15*d) + (11*b*cos(c + d*x)^2*sin(c + d*x))/(15*d) - (b*cos(c + d*x)
^4*sin(c + d*x))/(5*d)